On the role of chemical synapses in coupled neurons with noise 
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We examine the behavior in the presence of noise of an array of Morris-Lecar neurons coupled 
via chemical synapses. Special attention is devoted to comparing this behavior with the better 
known case of electrical coupling arising via gap junctions. In particular, our numerical simulations 
show that chemical synapses are more efficient than gap junctions in enhancing coherence at an 
optimal noise (what is known as array-enhanced coherence resonance): in the case of (nonlinear) 
chemical coupling, we observe a substantial increase in the stochastic coherence of the system, in 
comparison with (linear) electrical coupling. We interpret this qualitative difference between both 
types of coupling as arising from the fact that chemical synapses only act while the presynaptic 
neuron is spiking, whereas gap junctions connect the voltage of the two neurons at all times. This 
leads in the electrical coupling case to larger correlations during interspike time intervals which 
are detrimental to the array-enhanced coherence effect. Finally, we report on the existence of a 
system-size coherence resonance in this locally coupled system, exhibited by the average membrane 
potential of the array. 



are more efficient at enhancing coherence than gap junc- 
tions. We provide a qualitative explanation for this fact, 
paying particular attention to the roles of synchroniza- 
tion and of the correlation between neural dynamics in 
the time lapse between spikes. To that end, we exam- 
ine the effect of a linear pulsed coupling designed ad hoc 
for this purpose, which exhibits the optimal features of 
chemical coupling while still acting linearly on the mem- 
brane potential during spiking. 



I. INTRODUCTION 

Neurons are excitable devices that respond in a spiky 
manner to extrinsec stimuli. These stimuli can be pro- 
vided by external excitation, by noise, or by neighboring 
neurons in an extended system yj. In the absence of 
deterministic external driving, isolated neurons exhibit a 
spiking behavior purely induced by noise, with the pe- 
culiarity that the temporal coherence of the system in- 
creases for increasing noise up to a certain noise level, be- 
yond which coherence decreases again. Thus, an optimal 
amount of noise exists for which coherence is maximal. 
This phenomenon, which we call stochastic coherence (to 
stress the analogy with the better known stochastic reso- 
nance), is known in the literature as coherence resonance 
or internal stochastic resonance 0, Q • 

Recent studies have shown that, in extended arrays 
of neurons, coupling noticeably enhances the stochastic 
coherence effect 0, HI- This array- enhanced stochastic 
coherence (AESC) has been reported so far, up to our 
knowledge, only in the case of linear (diffusive) electrical 
coupling, mediated by gap junctions between the neu- 
rons 0. But another very important means of signal 
transmission between neurons is via chemical synapses, 
which provide a nonlinear pulsed coupling only when the 
presynaptic neuron is excited. It is thus of interest to ex- 
amine the effect of this kind of nonlinear coupling on the 
stochastic resonance effect described above. Our numer- 
ical results, detailed below, show that chemical synapses 



II. MODEL DESCRIPTION 
A. Neuron model 

We consider a one-dimensional array of neurons whose 
dynamical behavior is described by the Morris-Lecar 
model Q, 



dt 
dWj 
dt 



1 

Cn; 



(i) 

(2) 



where i — 1, . . . N index the neurons, and V* and Wi rep- 
resent the membrane potential and the fraction of open 
potasium channels, respectively. C m is the membrane 
capacitance per unit area, 7? pp is the external applied 
current, J- yn is the synaptic current, and the ionic cur- 
rent is given by 
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where g a (a = Ca, K, L) are the conductances and V® the 
resting potentials of the calcium, potassium and leaking 
channels, respectively. <f> is the decay rate of Wi, and we 
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define the following functions of the membrane potential: 

'V-Vmx ' 



neuron i is given by 



Moo(V) = \ 

Woo(V) = \ 

A(V) = cosh 



1 + tanh 

1 + tanh 
V-Vwi 



Vm2 

V-Vwi 

Vw2 



(4) 
(5) 
(6) 



where Vmi, Vm2, Vwi an d Vwi are constants to be spec- 
ified later. The last term in Eq. is a white Gaussian 
noise term of zero mean and amplitude Di, uncorrelated 
between different neurons. 

In the absence of noise, an isolated Morris-Lecar neu- 
ron shows a bifurcation to a limit cycle for increasing 
applied current 7 app Q • This bifurcation can be a saddle- 
node (type I) or a subcritical Hopf (type II) bifurcation, 
depending on the parameters. We chose this last option 
for the numerical calculations presented in this paper. 
The specific values of the parameters used are shown in 
tabled 9]. The equations were integrated using the Heun 
method [l(|, which is a second order Runge-Kutta algo- 
rithm for stochastic equations. 



B. Coupling scenarios 



Most of the studies done so far in the field of stochas- 
tic neural dynamics consider linear electrical coup ling 
thr ough g ap j unctions between adjacent neurons [3, HH, 
[3 EJ EJ EH • But the most common way used by neu- 
rons to transmit information is by means of nonlinear 
pulsed coupling through chemical synapses. In the fol- 
lowing paragraphs we discuss the modeling of the two 
types of coupling, both of which will be analyzed later. 



1. Linear diffusive coupling: gap junctions 

In this kind of coupling, the synaptic current is pro- 
portional to the membrane potential difference between 
a neuron and its neighbors, 



E 

jSncigh(i) 



(7) 



where Vj stands for the membrane potential of neuron i, 
the sum runs over the neighbors that feed that neuron, 
and g^ n is the conductance of the synaptic channel. 



2. Nonlinear pulsed coupling: chemical synapses 

In order of take into account the chemical nature of 
the synapses, we use the model proposed in 0] to couple 
the neurons. In this model, the synaptic current through 



E 

jencigh(j) 



(8) 



where the sum runs over the neighbors that feed neuron 
i, g^ yn is the conductance of the synaptic channel, Tj rep- 
resents the fraction of bound receptors, Vi is the postsy- 
naptic membrane potential, and E s is a parameter whose 
value determines the type of synapse (if larger than the 
rest potential, e.g. E s — mV, the synapse is excitatory; 
if smaller, e.g. E s = —80 mV, it is inhibitory). 

The fraction of bound receptors, rj, follows the equa- 
tion 



(9) 



where [T] 3 = 6»(Tj + r syn - t)6(t - Tj) is the concentra- 
tion of neurotransmitter released into the synaptic cleft, 
a and (3 are rise and decay time constants, respectively, 
and Tq is the time at which the presynaptic neuron j fires, 
which happens whenever the presynaptic membrane po- 
tential exceeds a predetermined threshold value, in our 
case chosen to be 10 mV. This thresholding mechanism 
lies at the origin of the nonlinear character of the chem- 
ical synaptic coupling, which contrasts with the linear 
nature of the diffusive electrical coupling of Eq. Q. The 
time during which the synaptic connection is active is 
given by T syjj . The values of the coupling parameters 
that we use |lg are specified in Table [I] . 



Parameter 


Value 


Cm 


5AiF/cm 2 


9K 


8/iS/cm 2 


9L 


2 fiS/cm 


gca 


4.4^S/cm 2 


Vk 


-80 mV 


Vl 


-60 mV 


Vca 


120 mV 


Vmi 


-1.2 mV 


Vmi 


18 mV 


Vwi 


2mV 


V\V2 


30 mV 




l/25s _1 


a 


0.5 ms^mivr 1 


P 


0.1ms -1 




(specified in each case) 


Tsyn 


(specified in each case) 


E s 


OmV 



TABLE I: Parameters used in this work. 
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III. THE BEHAVIOR OF N COUPLED 
NEURONS: CHEMICAL SYNAPSES VS GAP 
JUNCTIONS 

To obtain a first glimpse of the different effects of chem- 
ical and electrical synapses in neural dynamics, we begin 
by studying the behavior of two coupled neurons under 
the influence of increasing noise intensity D, for the two 
types of coupling. We fix the external current for the two 
neurons at 7 a PP = 46 mA, in such a way that they do not 
fire in absence of noise. The conductance in the synaptic 
channels (a measure of the coupling strength) was cho- 
sen so that the neurons fire synchronously at least for low 
noise intensities {g syn = 4 nS for chemical synapses and 
g syn = 1 nS for electrical coupling). We use independent 
realizations of the noise in each neuron but with the same 
amplitude D. 

In the deterministic case (D — 0) the system is in a 
quiescent state. At small but nonzero noise levels the 
neurons spike sparsely (albeit synchronously), producing 
an irregular sequence of well-separated spikes, as can be 
seen in the top row of Fig. ^ For an intermediate noise 
amplitude the spiking rate increases, as does the regular- 
ity of the spiking events (middle row in Fig.^l. Finally, 
for strong noise the spiking of the neurons becomes irreg- 
ular again (bottom row in Fig. ^| . The situation is quali- 
tatively the same for both types of coupling, even though 
when comparing the most optimal situation in each case, 
the coherence in the case of chemical coupling is larger 
than for electrical coupling (as we quantify below). We 
note that even though the two neurons fire synchronously 
in both cases, in the chemical-coupling case there exists 
a slight delay between spikes, and the inter-spike dynam- 
ics are different from each other (cf the solid and dashed 
lines in all plots of Fig.^l. In the electrical-coupling case, 
on the other hand, the time traces are basically identical 
at all times. This is our first indication that gap junc- 
tions are much more efficient than chemical synapses in 
leading to synchronization. 

As usual in neurophysiology, in order to quantify the 
behavior shown above we evaluate the time interval be- 
tween consecutive spikes, T p , as the main observable in 
our numerical simulations. In particular, we analyze the 
first two statistical moments of the distribution of T p , 
namely its mean value (T p ) and its normalized stan- 
dard deviation (also known as coefficient of variation) 
Rp = <Tp/ (T p ). 

Figure |21 plots these two quantities vs. noise strength 
for both chemical and electrical coupling in linear arrays 
of N = 1,2,10 and 30 neurons, coupled bidirectionally 
and with periodic boundary conditions. While for both 
types of coupling the average interspike interval (T p ) de- 
cays quickly with noise and levels off independently of 
the system size, the coefficient of variation shows clear 
differences between chemical and electrical coupling. In 
the former case, R p shows a clear trend as N varies, with 
both its minimum value and the corresponding optimal 
noise level decreasing steadily (and strongly) with in- 



chemical coupling 



electrical coupling 




time (ms) 



time (ms) 



FIG. 1: Membrane potential time traces of two neurons 
coupled via chemical synapses (left panels) and gap junc- 
tions (right panels) for three differents noise amplitudes: in 
the chemical-coupling case D — 0.5 mV/ms (top), D — 
3.5 mV/ms (middle) and D = 12 mV/ms (bottom); in the 
electrical-coupling case D — 1 mV/ms (top), D — 4 mV/ms 
(middle) and D = 12 mV/ms (bottom). The time series of 
both neurons are shown as solid and dashed lines (almost in- 
distinguishable in the top right panel). J app = 46 mA, and 
Tsyn = 1.5 ms in all cases. For chemical coupling, g syn — 4 nS 
(left), whereas for the electrical coupling g syn = 1 nS (right). 



chemical coupling 



electrical coupling 




D (mV/ms) 



D (mV/ms) 



FIG. 2: Mean time interval between spikes (T p ) (upper pan- 
els) and coefficient of variation R p = a v j (T p ) (lower pan- 
els) for the membrane potential of one neuron coupled via 
chemical synapses (left panels) and gap junctions (right pan- 
els). I app — 46 mA and r syn = 1.5 ms in all cases, whereas 
g syn = 4 nS for the chemical coupling (left), and g syn = 1 nS 
for the electrical coupling (right). 



creasing N. For electrical coupling, on the other case, 
the decrease is less pronounced. For a given system size, 
the coherence of the chemically coupled array is much 
better than that of the electrically coupled array. 

It could be argued that the difference in behavior be- 
tween the two types of coupling is due to the different 
coupling coefficients g syn used. But an analysis of the 
influence of this parameter in each model confirms that 
chemical coupling is overall more efficient than the elec- 
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FIG. 3: Minimum coefficient of variation R p — a p /(T p ) of the 
local (top) and average (bottom) membrane potential for two 
coupled neurons versus coupling strength, for both chemical 
and electrical coupling. I apP = 46 mA in all cases. 

trical one in enhancing coherence. This is evidenced by 
Fig. Ola) , which plots the minimum coefficient of varia- 
tion R p (i.e. its value for optimal noise) versus the cou- 
pling coefficient g syn for the two types of coupling in the 
2-neuron case. It is clear that optimal coherence under 
electrical coupling is only attainable for a narrow range 
of coupling strengths. For small and large coupling lev- 
els coherence is low, and only for a small intermediate 
range of coupling strengths coherence is enhanced. On 
the other hand, for chemical coupling the range of cou- 
pling strengths for which coherence is enhanced is much 
larger, apparently not being bounded from above. 



IV. SYNCHRONIZATION AND ARRAY 
ENHANCEMENT 

When choosing the right values of the coupling coeffi- 
cient that make chemical and electrical coupling compa- 
rable, one also needs to take into account the synchro- 
nization between the different elements in the array. This 
can be accomplished via the average membrane potential 
Vkvg over all the array. In the case of perfect synchroniza- 
tion, this quantity will exhibit spikes identical to, and 
simultaneous with, those of the individual neurons. In 
the case of partial synchronization, spikes in the average 
potential (defined beyond a given threshold, in our case 
10 mV) only occur when a sufficient number of neurons 
fire within a time window smaller than the spike width. 
Given these considerations, we determine the interspike 
interval series of the average potential and compute its 
statistical properties, namely its mean (T p ) and its nor- 
malized standard deviation R p = a p /{T p ). 

Figure|3{b) plots the minimum value of R p (for optimal 
noise) corresponding to the average potential, for increas- 
ing coupling strengths. On the basis of these results, we 



can establish that the best choices for coupling strengths 
are <? syn = 1 nS for electrical coupling and <? syn > 4 nS 
for chemical. These values are valid not only for two 
neurons, but for all analyzed system sizes. 

We have seen so far that chemical coupling is more 
effective at enhancing stochastic coherence than electri- 
cal coupling. We now discuss the relationship between 
array-enhanced coherence and synchronization. Array- 
enhanced coherence is absent for low levels of synchro- 
nization, because in that case every neuron behaves as 
basically uncoupled from all others, so that no construc- 
tive effects of coupling (such as array-enhanced coher- 
ence) can arise. In the opposite limit, large synchroniza- 
tion levels are not useful in promoting coherence either, 
since in that case the whole system behaves as a single 
neuron. However, in the intermediate regime of imper- 
fect synchronization, neighboring neurons can "remind" 
each other to spike at the right times (i.e., right after the 
refractory period) for an optimal noise level, so that co- 
herence is globally enhanced. This cannot happen if the 
neurons are perfectly synchronized. 

We now test the efficiency of the two coupling schemes 
in leading to synchronization. We can already expect, 
from the functional form of the corresponding coupling 
term, Eqs. that chemical synapses will never lead 

to perfect isochronous synchronization in neuronal ar- 
rays. To begin with, the dynamics of the fraction Tj of 
bound receptors introduces a delay in the interneuronal 
communication that is absent in the electrical coupling 
case, Eq. {7J. The diffusive form of the latter, further- 
more, is compatible with the identical synchronization 
solution, Vi(t) = Vj(t) Vi,j, whereas this is not so for 
chemical coupling. 

To check the previous hypothesis, Fig. 0] compares the 
mean interspike interval and coefficient of variation of the 
average potential for the two types of coupling. This fig- 
ure should be contrasted with the corresponding results 
for the local membrane potential, Fig. [21 Again, the situ- 
ation is qualitatively similar in both cases, although there 
are quantitative differences. The mean time interval be- 
tween spikes of the average potential increases noticeably 
for large system sizes (c/ stars in the top row of Fig. 0}, 
but much more for chemical coupling. Also, the coeffi- 
cient of variation of the interspike interval series becomes 
more degraded (i.e. the system response is worst at the 
same system sizes) in the chemical synapse case. These 
facts indicate that synchronization is worse for chemical 
than for electrical coupling. 

In order to confirm the previous conclusions, we follow 
how spikes propagate through the array by means of the 
spike diagram shown in Fig. [S] In this diagram, we plot 
vertical lines every time that one neuron spikes. In this 
kind of diagram synchronized spiking is clearly identified 
by vertical alignment of the spike markers. 

Figure O plots the spike diagram for the two types of 
coupling discussed so far and for N = 2, 5, 10 and 30 
neurons. By comparing the two columns we can observe 
that the system with chemical synapses has a charateris- 
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chemical coupling electrical coupling 




D (mV/ms) 



FIG. 4: Mean time between spikes (T p ) (upper panels) and 
coefficient of variation R p = a v j {T p } (lower panels) for the 
average membrane potential and increasing system sizes. Pa- 
rameters are those of Fig. 
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FIG. 5: Spike diagram for TV neurons coupled via chemical 
synapses (g syn — 4 nS left panels) and gap junctions (<? syn = 
1 nS rigth panels) for 2, 5, 10 and 30 neurons. I apP = 46 mA 
and D = 3 mV/ms in all cases. 



tic delay time in the propagation of spikes between neu- 
rons larger than the corresponding to the electrical cou- 
pling. Due to this delay, synchronization is rapidly lost 
in the case of chemical coupling. For electrical coupling 
the spikes also show slight propagation delays, but much 
smaller than in the chemical case. With the values of 
coupling constants chosen here, both systems lose syn- 
chronicity at similar sizes. If we increased the value of 
the g syn in the electrical coupling case, the system would 
remain synchronized for higher system sizes. This would 
not happen in the chemical case, because there the delay 
is intrinsic to the synapses and does not depend on the 
coupling strength. This is a fundamental, qualitative dif- 
ference between chemical and electrical coupling, which 
is worth being highlighted. 



V. NONLINEAR VERSUS PULSED COUPLING 

There are two main differences between coupling via 
chemical synapses and gap junctions. First, electrical 
coupling occurs continuously whereas chemical synapses 
are only active when the presynaptic neuron spikes. Sec- 
ond, chemical coupling is intrinsically nonlinear, because 
beyond a certain threshold for the presynaptic neuron, 
the coupling signal has always the same shape. Which 
one of these two features leads to enhancement of stochas- 
tic coherence reported above for the chemical coupling? 
We postulate that the relevant feature that makes the 
chemical synapses more efficient to help neurons to fire 
is that chemical coupling is only effective when the presy- 
naptic neuron fires. Otherwise, two neighboring neurons 
are uncoupled. 

To test this hypothesis we couple the neurons with an 
artificial synaptic current proportional to the difference 
between the voltages of the coupled membranes (as in 
the usual linear coupling described in Sect. IIIBl"|) . but 
just during a time r syn after the presynaptic neuron has 
fired. This model, that we can call linear pulsed coupling, 
is a hybrid of the two models used above: the coupling is 
linear with the membrane potential, but the neurons are 
uncoupled as far as the presynaptic neuron does not fire. 

The synaptic current affecting neuron i due to the in- 
teraction with its neighboors can be written as 

7f n = ]T 9(Tj+ v -t)9(t-T^rr j (V i -^), 

jeneigh(i) 

(10) 

where 9{T^ + r syn -t)6{t- T J ) is a pulse of width T syn 
that turns on when the membrane potential of neuron j 
is larger that a certain threshold (lOmV in the present 
case). With this model, a neuron is uncoupled from its 
neighbors if they are silent, and coupled during a time 
T syn after one of them fires. 

We also calculate in this case the average time between 
consecutive spikes, (T p ) and its normalized standard de- 
viation R p — o p /(T p ), both for the spikes produced by 
the membrane potential of one neuron and for the av- 
erage membrane potential. The results are plotted in 
Fig. El 

Analyzing the results for the local membrane potential 
(left panel of Fig.^J, we can see that the behavior in this 
case is similar to (even better than) the case of full chem- 
ical coupling: the response of a neuron to a purely noisy 
excitation is clearly enhanced already when only two neu- 
rons are coupled, and rapidly improves with system size, 
saturating after N — 10 neurons. Also, the coefficient 
of variation R p is as low as for chemical synapses. This 
happens even though coupling is linear when active. 

If we analyze the behavior of the average potential 
(right panel of Fig. EJ, we can see that synchronization 
holds at even larger system sizes than for the previous 
natural coupling schemes. In this hybrid model, the 
(Vi — Vj) term enhances synchronization for large cou- 
pling strengths, but without losing enhancement as in 
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local membrane potential 



average membrane potential 




D (mV/ms) 



D (mV/ms) 



FIG. 6: Mean time between spikes (T p ) (upper plot) and co- 
efficient of variation R p = ct p /{T p ) (lower plot) for the mem- 
brane potential of one neuron (lef panel) and for the average 
membrane potential (right panel). The neurons are coupled 
linearly according with expression 1101 Parameters are those 
of Fig. El (&,,,» =4nS). 




FIG. 7: Coefficient of variation R p for the membrane poten- 
tial of one neuron as a function of system size for chemical 
(D — 1 mV/ms) and electrical (D — 3 mV/ms) coupling. 
gayn _ ^ n g f Qr chgjjjigai coupling and g syn = 1 nS for the 
electrical coupling. 



VII. DISCUSSION 



the linear coupling model, due to the inactivation of the 
coupling in the time elapsed between spikes. 

In summary, array-enhanced stochastic coherence is 
present for linear pulsed coupling. The fact that, in this 
hybrid type of coupling, this term only turns on when 
presynaptic neurons fire, supports the hypothesis that 
the most important feature underneath the efficiency of 
the coupling is that the neurons remains uncoupled in 
absence of firing. 



VI. SYSTEM-SIZE COHERENCE RESONANCE 

The behavior exhibited in Fig. 0] by the coefficient of 
variation R p of the average membrane potential in the 
chemical-coupling case displays a relevant feature: for 
small enough noise levels (such as D = 1 mV/ms), as 
the system size increases the regularity of the system in- 
creases (i.e. the coefficient of variation decreases), only 
to degrade again for even larger system sizes. This is a 
system-size resonance effect. System-size coherence res- 
onance has been recently reported in globally coupled 
excitable elements and in arravs of neurons interact- 
ing via linear diffusive coupling jl8| . Here we report it 
for chemically coupled neurons. 

Figure [7| compares the effect of system size on the 
spike-train coherence for chemically and electrically cou- 
pled neurons. Whereas in both cases, the lowest coher- 
ence is found for an optimal system size, the effect is 
more pronounced, and occurs at the same system size for 
the values of g syn chosen. 



Recent studies have shown that gap junctions are more 
efficient than chemical synapses in leading to synchro- 
nization . The evidence presented in this paper shows 
that chemical synapses are substantially more efficient 
than gap junctions in enhancing stochastic coherence in 
coupled neuron systems. This difference in efficiency 
stems from the fact that electrical coupling enhances cor- 
relations among neurons during interspike time intervals, 
which prevents neighboring elements to "remind" each 
other to fire at the right time. Chemical coupling, on the 
other hand, has an intrinsic transmission delay and an in- 
trinsic electrical isolation that degrades this correlation, 
making that cooperative effect possible. This is also re- 
flected in the different behavior of both types of coupling 
versus coupling strength (Fig. |3J) : for electrical coupling 
coherence is low only for a narrow window of coupling 
strengths; for chemical coupling, on the other hand, low 
coherence arises even if coupling is large. Clearly this 
difference arises from the fact that for diffusive coupling, 
large coupling levels lead to high synchronization; for 
chemical coupling, on the other hand, synchronization 
does not improve even if coupling strength increases, due 
to the intrinsic delay. 

Interneuron communication via chemical synapses is 
ubiquitous. Its role must therefore be properly assessed 
when studying stochastic effects in neuronal dynamics. 
The results presented here show that this type of coupling 
is in fact beneficial for the system coherence. 
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